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Abstract. We present the detailed spectroscopic analysis of 72 evolved stars, which were previously studied for accurate 
radial velocity variations. Using one Hyades giant and another well studied star as the reference abundance, we determine the 
[Fe/H] for the whole sample. These metallicities, together with the T e g values and the absolute V-band magnitude derived from 
Hipparcos parallaxes, are used to estimate basic stellar parameters (ages, masses, radii, (B-V)q and logg) using theoretical 
isochrones and a Bayesian estimation method. The (B-V)o values so estimated turn out to be in excellent agreement (to 
within ~ 0.05 mag) with the observed (B-V), confirming the reliability of the r c ff-(B-V)o relation used in the isochrones. 
On the other hand, the estimated logg values are typically 0.2 dex lower than those derived from spectroscopy; this effect has 
a negligible impact on [Fe/H] determinations. The estimated diameters 9 have been compared with limb darkening-corrected 
ones measured with independent methods, finding an agreement better than 0.3 mas within the 1 < 6 < 10 mas interval (or, 
alternatively, finding mean differences of just 6 %). We derive the age-metallicity relation for the solar neighborhood; for the 
first time to our knowledge, such a relation has been derived from observations of field giants rather than from open clusters 
and field dwarfs and subdwarfs. The age-metallicity relation is characterized by close-to-solar metallicities for stars younger 
than ~ 4 Gyr, and by a large [Fe/H] spread with a trend towards lower metallicities for higher ages. In disagreement with other 
studies, we find that the [Fe/H] dispersion of young stars (less than 1 Gyr) is comparable to the observational errors, indicating 
that stars in the solar neighbourhood are formed from interstellar matter of quite homogeneous chemical composition. The 
three giants of our sample which have been proposed to host planets are not metal rich; this result is at odds with those for 
main sequence stars. However, two of these stars have masses much larger than a solar mass so we may be sampling a different 
stellar population from most radial velocity searches for extrasolar planets. We also confirm the previous indication that the 
radial velocity variability tends to increase along the RGB, and in particular with the stellar radius. 



1. Introduction 

It has recently become evident from high accuracy radial veloc- 
ity (RV) measurements that late-type (G and K) giant stars are 
RV variables (see e.g. Hatzes & Cochran 1993, 1994; Setiawan 
et al. 2003, 2004). These variations occur on two greatly differ- 
ent timescales: short term variability with periods in the range 
2-10 days, and long term variations with periods greater than 
several hundreds of days. The short term variations are due to 
stellar oscillations. The cause of the long term variability is not 
clear and at least three mechanisms, i.e., low mass compan- 
ions, pulsations and surface activity, have been proposed. The 
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first results of our long term study of the nature of the vari- 
ability of K giants have shown that all three mechanisms are 
likely contributors to long term RV variability, although their 
dependence on the star's fundamental characteristics remain 
unknown (Setiawan et al. 2003, 2004). 

The purpose of this paper is to accurately determine the 
radii, temperatures, masses, and chemical composition of the 
stars of our sample, in order to understand better how RV vari- 
ability and stellar characteristics are related. We are particularly 
interested in determining to what extent the stars found to host 
planetary system class bodies (Setiawan et al. 2003a, 2005) ex- 
hibit high metallicity, as has been found for the dwarfs hosting 
giant exoplanets ( Gonzales 1997; Santos et al. 2004). 

The present sample is composed of the stars published by 
Setiawan et al. (2004). The selection criteria and the observa- 
tions were explained there. We recall that observations were 
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obtained with the FEROS spectrograph at the ESO 1 .5m tele- 
scope (Kaufer et al. 1999), with a resolving power of 50000 
and a signal-to-noise ratio exceeding 150 in the red part of the 
spectra. 

2. Atmospheric parameters determination 

The spectroscopic analysis has been made in LTE, using a mod- 
ified version of Spite's (1967) code. MARCS plane parallel 
atmosphere models are used. Gustafsson et al. (1975) models 
were used for the giants, while Edvardsson et al. (1993) models 
were used for the few dwarfs of our sample. No major spurious 
effects are expected by using those two sets of model atmo- 
spheres (Pasquini et al. 2004). Equivalent widths of spectral 
lines were measured using the DAOSPEC package (Pancino & 
Stetson 2005). The spectra of several stars were cross checked 
by measuring the equivalent widths with MIDAS, resulting in 
a very good agreement of the two methods. The line list and 
corresponding atomic data were those adopted by Pasquini et 
al. (2004). 

Atmospheric parameters (effective temperatures, surface 
gravities, microturbulence velocities, and metallicities) have 
been obtained using an iterative and self-consistent procedure. 
The atmospheric parameters were determined from the spectro- 
scopic data in the conventional way. The effective temperature, 
r e fF, is determined by imposing that the Fe I abundance does 
not depend on the excitation potentials of the lines. The mi- 
croturbulence velocity is determined by imposing that the Fe I 
abundance is independent of the line equivalent widths. The Fe 
I/Fe II ionization equilibrium has been used to determine the 
surface gravity g. The method used here is described in detail 
by d el Peloso et al. (2005)| with the difference that temperature 
and gravity are also free parameters in the present work. The 
initial values for T e g and log g to start the process were derived 
in two steps. The first step involves the (B-V) photometry and 
parallax from Hipparcos, assuming for all stars solar metallic - 
ity and a 1.5 M mass, and adopting the |Alonso et al. (1 999) 
scale. The initial microturbulence velocity was 1.0 kms -1 . In 
the second step, the values for logg were obtained from the 
Girardi et al. (2000) evolutionary tracks by interpolation, using 
the previously determined values for T e g, abundance and My ■ 

2.1. Effective temperature 

The {V-K) index, which is used by many authors, is proba- 
bly the best photometric T e g indicator for G-K giants (Plez 
et al. 1992, Ramirez & Melendez 2004). The only source of 
/f-band photometry for our sample is the 2MASS Catalog 
|(Cutri et al 2003)| but its authors caution that stars with < 3 
are saturated, and the error of their colors is larger. Many 
stars of our sample belong to this class and the rest is not 
much fainter. We compared the r e ff values obtained from 
(V-K) with those from (B-V) to verify whether we can 
use the 2MASS catalog colors to determine the effective 
temperature of our sample. For both indices we used the 
|Alonso et al. (199 9)1 calibrations for the stars with \ogg < 
3.5 and the |Alonso et al. (1996)| calibration for the remain- 
ing dwarfs. The 2MASS K s magnitude was converted to 



the CTS system ( |Alonso etal. (1998)) via the CIT system 
|(Cutri et al 2003)| The (V-K) indices of the stars with log g > 
3.5 were converted to the Johnson system using the relation 
given by |Alonso et al. (1998)| The [Fe/H] values needed to ap- 
ply those relations were taken from Tabled 

The results are presented Fig. \l\ It is evident that the 
TeSty-K) vs. r e ff (B _ V ) relation has a much larger dispersion than 
Trffspec vs - T e fi(B-v)- The dispersion of the relation using J-H 
and J-K is greater. These findings confirm that the colors of 
2MASS are unsuitable to determine the r e ff of bright stars, in 
particular of our sample. On the other hand, we see in Fig. ^ 
that r e ff sp ec is m good agreement with T e s(B-v) f° r stars with 
4200 < r e ff < 5200, while there is less agreement for cooler 
and hotter stars. This result is expected for cold stars - because 
is known that (B-V) is not a good T e s indicator for them - 
but it is unexpected for those hotter stars. Noting that the hotter 
stars of our sample are dwarfs, we compared the r e (j which we 
determined from spectroscopy and from the (B-V) index with 
the values given in |del Peloso et a L(2005)| for a sample o f G- 
K dwarfs. The temperatures of the ldelPefoso et al. (2 005) stars 
were carefully determined using criteria different from the Fe 
I excitation equilibrium. The agreement between their and our 
results is better for our spectroscopy values than for the (B-V) 
values. Note that other authors also found differences between 
spectroscopic and photometric 7*^ for dwarfs of the same or- 
der or larger than the differences we found (e.g. Reddy et al. 
2003; Santos et al. 2004; Ramirez & Melendez 2004; Luck & 
Heiter 2005). Moreover, for T e fi < 4000A" there is a saturation 
of (B-V) and this colour index is no longer useful (see Alonso 
et al. 1999 for a discussion about this point). For these reasons, 
we will adopt our spectroscopic temperatures in the followin 
analysis. 

2.2. Microturbulence and iron abundance 

The determination of microturbulence in giants may be diffi- 
cult, as discussed in detail in previous papers (e.g. McWilliam 
1990). In order to test the sensitivity of our analysis with regard 
to this parameter, we have performed the process described 
above by starting with two very different values of microtur- 
bulence (e.g. 1 and 2.5 kms -1 ) and letting the system converge 
freely. We found for almost all stars a very good convergence, 
resulting in the same final values for all parameters, including 
microturbulence. Only for the four coolest stars of the sample 
does the analysis converge to microturbulence velocities that 
differ by up to 0.2 km s _1 , and to Fe abundances that differ by 
up to 0. 1 dex. These results were obtained using the same line 
lists and equivalent widths. This result indicates that the so- 
lution space may be degenerate for the coolest stars, and that 
fairly large errors may intrinsically be present. We also notice 
that for these stars, we may be exceeding the validity limit of 
the grid of our atmosphere models. The atmospheric parame- 
ters (effective temperature, gravity, metallicity and microturbu- 
lence) which we determined from the spectral analysis are pre- 
sented in Tabled In a coming paper (da Silva et al., in prepara- 
tion) we will present and analyze the abundances of some other 
elements and we will publish the equivalent widths and atomic 
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Table 1. Retrieved atmospheric parameters for the entire star sample. Note that metallicities are 0.07 dex larger in comparison 
to the scale derived independently from the analysis of HD 2737 1 and HD 1 1 3226; see text. Column 1 : HD number; Column 2: 
Spectroscopic effective temperature; Column 3: Iron content (normalized to the Sun); Column 4: Spectroscopic gravity; Column 
5: Microturbulence (kms -1 ). 
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data (log gf and excitation potential) of all the lines used in our 
analysis, for all sample stars. 

We still have to determine the zero point of the metallicity 
scale. Pasquini et al. (2004) used a normalization to the solar 
spectrum, because their sample contains many solar-type stars. 
Since our sample is mainly composed of evolved stars, we pre- 
fer to fix the abundance scale by using giants that have been 
well studied in the literature. Two stars of our sample are par- 
ticularly well studied, with several high quality entries in the 
Cayrel de Strobel et al. (2001) catalogue,namely HD 113226 
(e Vir) and HD 27371 (y Tau), a Hyades giant (the literature 
data for another Hyades giant in our sample, HD 27697, are 
much less constrained). 

Averaging all entries of the Cayrel de Strobel et al. (2001) 
catalogue forHD 113226, we obtain [Fe/H] = 0.11 and r e ff = 
5048 K, while for HD 27371 we obtain [Fe/H] =0.11 and 
T e g = 4967 K (the two values of Komaro in the catalogue 
for HD 27371 were not considered). Our results from Tabled 
are [Fe/H] =0.16 and r eff = 5086 K, and [Fe/H] = 0.20 
and T e ff = 5030 K for the two stars. We conclude that it is 
very likely that the metallicity scale of Table ^ is too high by 
0.07 dex in [Fe/H], and we will apply this correction to all our 
data in the following. This correction should therefore be ap- 
plied to all the entries of Table[0to derive the correct metallic- 
ity. Table [6] shows the corrected values. We note also that our 
spectroscopic temperatures for HD 113226 and HD 27371 are 
slightly higher (about 50 K) than the average values reported in 



the literature. Both results are in very good agreement with the 
results of Pasquini et al. (2004), who found [Fe/H] = 0.06 for 
the analysis of the UVES solar spectrum, and who derived sys- 
tematic differences between the spectroscopic and photometric 
temperatures of IC 4651 stars. Note also that the largest dif- 
ferences between the average and individual [Fe/H] values in 
the catalogue are 0.10 dex for HD 113226 and 0.09 dex for 
HD 27371, both of which are larger than the discrepancy with 
the present analysis. 

We notice also that another star in our sample has more 
than two entries in the Cayrel de Strobel et al. (2001) cata- 
logue: the moderately metal poor star HD 18907, which is at 
the low metalliciy end of our sample. The average literature 
value for this star is [Fe/H] = -0.69, while our zero-point cor- 
rected value is [Fe/H] = -0.54. Although we do not expect a 
strong dependence of the zero point on the metallicity, if we 
were to include this star in our set of calibrators, we should 
apply a correction of 0.10 dex instead of 0.07. In our analysis, 
HD 18907 has a very low microturbulence (0.9 kms~'), which 
is lower than for any other star, providing an explanation of the 
discrepancy with the literature. Although a second analysis of 
this star did not reveal any suspicious effect, we preferred not 
to use it to determine our zero-point. 

We could redraw Fig. 1 with the T e gfB-v) an d ^efftv-AO calcu- 
lated using the corrected values of [Fe/H], given in Table[6] but 
the r e ff (B _y) calibration is not very sensitive to the metallicity, 
and T e ff (V _ K) is even less so. For the stars cooler then 4500 K 
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Fig. 1. r e ff values obtained from (B—V) compared with the val- 
ues obtained from the spectroscopic analyses (left panel), and 
compared with r e ff obtained from (V-K) (right panel). Giants 
are represented by circles and dwarfs by squares in both panels. 
The upper panels show the differences between the values 

there is no difference at all between the T e ff found from the two 
[Fe/H] values. The largest difference in our sample is 31 K, 
found for the hottest star HD 26923, from 7 , e ff( B _v ) , T e ff being 
lower for the correct (lower) metallicity. 

We compare the [Fe/H] and T e g values of our analysis with 
those from the Cayrel de Strobel et al. (2001) catalogue in 
Fig. |2 Stars with more than one entry in the catalogue are 
shown as empty squares. Most of the entries are from the work 
by Mc William (1990), and stars with only that entry are shown 
with filled squares. The results of the comparison are given in 
Table|5] and they show that the agreement with most data in the 
literature is quite good, and that we tend to systematically re- 
trieve higher abundances (by 0.1 dex) than McWilliam (1990). 
We somewhat expected such a result, because the method cho- 
sen by McWilliam to determine microturbulence tends to result 
in rather high values of However, his estimate of [Fe/H] for 
the Hyades giant HD 27371 is [Fe/H] = -0.02, which is much 
lower than our value, and lower than what is considered the 
best estimated abundance for this cluster. We would therefore 
expect a discrepancy between our and the McWilliam results 



Table 2. Results of the comparison of our stellar parameters 
with those of the Cayrel de Strobel et al. (2001) catalogue. All 
differences represent our values minus the literature values. The 
third and fourth rows present mean and standard deviation of 
the difference with respect to McWilliam (1990) only. 
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of about 0.13 dex, which is very close to our results. Note also 
that the values in Table |3 are not independent, since our av- 
erage computations from the literature include the McWilliam 
results. 

2.3. Iron abundance error estimate 

It is difficult to provide realistic error estimates for abundance 
and stellar parameter determinations, such as r e jf and gravity. 
We intend to be very careful in estimating these errors, because 
it is our goal to use the parameters to derive stellar masses and 
radii from evolutionary tracks. Any error in the fundamental 
(spectroscopic) parameters would make the determination of 
the derived stellar parameters more uncertain. 

The direct comparison of our results with those of other 
authors is shown in Tabled Our Fe content is on average sys- 
tematically too high by 0.07 dex, the temperature is too high 
by about 40 K and gravity is too high by about 0.13 dex. A 
systematic error of such a magnitude would not influence our 
measurements or conclusions in a significant way. However, 
the scatter about the mean discrepancies is in general more pro- 
nounced, and we believe that the scatter represents a more re- 
alistic, albeit somewhat pessimistic, estimate of the uncertainty 
of our retrieved parameters. 7^ shows a scatter of 70 K. The 
scatter of logg is up to 0.2 dex and quite large. The scatter of 
[Fe/H] is 0.1 dex. 

An estimate of the internal error in our [Fe/H] determi- 
nation, which is produced mainly by the uncertainties in the 
equivalent width measured and in the log gf is more impor- 
tant than the external error of ~ 0.1 dex (see next section). 
To estimate our internal errors in the determined metallicities, 
we examined two giants of our sample in more details, one 
among the coolest and the other among the hottest stars, namely 
HD 1 1 1884 and HD 36189. We changed each of the input pa- 
rameters (r e jf , log g and of our spectroscopic analysis in turn, 
using the values found above for the scatters (70 K, 0.2 dex and 
0. 1 dex). We proceeded in two ways: first, the other parameters 
were kept fixed; second, which is more in agreement with our 
method, the others parameters were left to adjust freely. The re- 
sulting changes of [Fe/H] were monitored and Table |3]presents 
the results. Note that step one is useful as a test because, in the 
method employed, when one parameter is changed, the others 
also change. An error in one parameter will be reflected in the 
other parameters. Usually the errors on the equivalent width 
measurements will produce greater dispersions in the diagrams 
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Fig. 2. Comparison of our results (abscissa) and literature values (ordinate) for [Fe/H], logg and T e g as retrieved from Cayrel 
de Strobel et al. (2001). Open squares represent stars with more than one determination in the literature; filled squares represent 
values from McWilliam (1990). 



Table 3. Sensitivity of [Fe/H] to changes of other stellar pa- Table 4. Sensitivity of stellar parameters to changes of equiva- 
rameters for two stars at the cool and hot ends of our sample. lent widths for two stars at the cool and hot ends of our sample. 
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3. Stellar parameters 



The bulk of our sample consists of giants and subgiants in 
the Hipparcos and Tycho catalogues (ESA 1997) with paral- 
laxes given with an accuracy better than 10% . This means 
that absolute My magnitudes are known with an accuracy of 
o"m v < 0.21 mag, whereas the apparent (B—V) colours (in the 
Johnson system) are also known quite precisely. The typical 
(B — V) error value given in the catalog is % 0.005 mag except 
for a few stars classified as variables. Moreover, reddening is 
expected to be negligible for this nearby sample, so that we 
could initially assume (B-V) = (B-V). Fig.|3]shows the sam- 
ple in the My vs. (B-V) diagram, with error bars included. 

Giant stars are well known to suffer the so-called age- 
metallicity degeneracy: old metal-poor stars occupy the same 
region of the color-absolute magnitude diagram (CMD) as 
young metal-rich objects. By having measured the metallicity 
of our sample stars, it should be possible to resolve this de- 
generacy and to estimate stellar ages from the position in the 
CMD. Some degeneracy will still remain, for instance we can- 



using them as a criterion, causing uncertainties in these quan- 
tities. The continuum placement, for instance, will be slightly 
higher in some regions and slightly lower in others. To illus- 
trate the influence of the equivalent width errors on our results, 
but stressing that this is not a real situation, we tested what 
happens when all those quantities are too large or too small by 
a constant factor. In our example we use 3%, the results are 
shown in Table |4] The errors are not symmetric and they are 
much enlarged for the hotter than the cooler stars. Fortunately, 
the continuum placement for the hotter stars are much more 
precise and a error of 3% is unlikely. Given these results, we 
will assume an internal error of 0.05 dex in our [Fe/H] deter- 
mination. 
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we compute the probability of the observed star to belong 
to this section 



B-V 



Fig. 3. Our sample in the color-absolute magnitude diagram. 
Error bars in My are derived from the parallax error of the 
Hipparcos catalog (<tm v = 2.17cr jr /;r), whereas errors in (B-V) 
are assumed to be the Hipparcos value or o~(b-v) = 0.005 mag, 
whichever is larger. The few stars with large (B-V) error 
bars are either known variables or - as for HD 124822, with 
0.3 mag-wide error bars - suspicious entries in the Hipparcos 
catalog. 



not easily distinguish between first-ascent RGB and post He- 
flash stars, or between RGB and early- AGB stars. 

3.1. Method for parameter estimation 

We implemented a method to derive the most likely intrinsic 
properties of a star by means of a comparison with a library 
of theoretical stellar isochrones, based on the ones by Girardi 
et al. (2000) and using the transformations to BV photometry 
described in Girardi et al. (2002). We adopt a slightly mod- 
ified version of the Bayesian estimation method idealized by 
J0rgensen & Lindegren (2005, see also Nordstrom et al. 2004), 
which is designed to avoid statistical biases and to take er- 
ror estimates of all observed quantities into consideration. Our 
goal is the derivation of complete probability distribution func- 
tions (PDF) separately for each stellar property x under study 
- where x can be for instance the stellar age t, mass M, surface 
gravity g, radius R, etc. The method works as follows: 

Given a star observed to have My + <tm v , T e ff ± cr rclr , and 

[Fe/H] ± CT[Fe/H], 

1 . we consider a small section of an isochrone of metallicity 
[Fe/H]' and age t' , corresponding to an interval of initial 
masses [Mi l ,M 2 ] and with mean properties M' T e s' and 
x'\ 



P 12 (x') 



J 



0(M;)dM; exp 



(log r eff - log T eS ') 2 



(M v - M v 'f 



cr 



(1) 



Mi 



log T e g 

where the first term represents the relative number of stars 
populating the [Mi l ,M 2 ] interval according to the initial 
mass function </>(M;), and the second term represents the 
probability that the observed My and T e ff correspond to the 
theoretical values, for the case that the observational errors 
have Gaussian distributions; 

3. we sum over P n (x') to obtain a cumulative histogram of 
P(x); 

4. we integrate over the entire isochrone; 

5. we loop over all possible [Fe/H]' values, now weighting the 
P l2 (x') terms with a Gaussian of mean [Fe/H] and standard 
deviation cr [Fe / H ]; 

6. we loop over all possible ?' values, weighting the P l2 (x') 
terms with a fiat distribution of ages; 

7. we plot the cumulative PDFs P(x) and compute their basic 
statistical parameters like mean, median, variance, etc. 

The oldest adopted stellar age is 12 Gyr, corresponding to an 
initial mass of about 0.9 M . Stellar masses as low as 0.7 M 
may be present in our PDFs, since mass loss along the RGB is 
taken into account in the Girardi et al. (2000) isochrones. 

The method implicitly assumes that theoretical models pro- 
vide a reliable description of the way stars of different mass, 
metallicity, and evolutionary stage distribute along the red gi- 
ant region of the CMD. This assumption is reasonable consid- 
ering the wide use - and consequent testing - of these models 
in the interpretation of star cluster data. Also, this assumption 
can be partially verified a posteriori, by means of a few checks 
discussed below. 

We take the logarithms of the age f, mass M, surface radius 
R and gravity g, together with the colour (B-V)q, as the pa- 
rameters x to be determined by our analysis. The reason to deal 
with logarithms is that their changes scale more or less linearly 
with changes in our basic observables - namely the absolute 
magnitude, log r e ff, and [Fe/H]. This choice of variables is ex- 
pected to lead, at least in the simplest cases, to almost symmet- 
ric and Gaussian-like PDFs. 

In comparison with Nordstrom et al.'s (2004) work, the 
mathematical formulation we adopt is essentially the same, ex- 
cept that 

- we use log T e s instead of T e tf in Eq.^ f° r the reason stated 
above, 

- we do not apply any correction for the a-enhancement of 
metal-poor stars, hence considering all stars in our sample 
to have scaled solar metallicities. In fact, few stars in our 
sample are expected to exhibit a-enhancement at any sig- 
nificant level. 

- we extend the method to the derivation of several stellar 
parameters, whereas Nordstrom et al. (2004) were limited 
to age and mass, 
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- we apply the method to red giants with spectroscopic 
metallicities, whereas Nordstrom et al. (2004) have a sam- 
ple composed mostly of main sequence stars with metallic- 
ities derived from uvby/3 photometry. 

Therefore, the empirical tests and insights provided by our 
sample will be different from, and hopefully complementary 
to, those already presented by other authors. 

We applied this method earlier, but assuming (B-V)o as 
the observable and T e ff as a parameter to be searched for. This 
means that T e g is replaced by (B—V) in Eq.[Q and vice-versa. 
Results of this early experiment are presented in Girardi et al. 
(2006), and are equivalent to the ones presented here with re- 
spect to the estimates of masses, ages, \ogg and radii. The use 
of a colour instead of the observed T e ff in Eq. [2 could be an 
interesting alternative for studies of red giants for which spec- 
troscopic r e jf is not accurate enough - e.g. because the avail- 
able spectra cover a too limited range and the [Fe/H] analysis 
is based on too small a number spectral lines. T e ff could then be 
derived from the photometry via a r e ff-color relation. We refer 
to Girardi et al. (2006) for a discussion of this point. 

The implementation of this method will soon be 
made available, via an interactive web form, at the URL 
http : //web . oapd. inaf . it/lgirardi 

3.2. Examples of PDFs 

Since one of the main achievements of this work is the applica- 
tion of the J0rgensen & Lindegren (2005) method to a sample 
of giants, it is important to illustrate the different situations we 
encounter for stars along the RGB. Figures |4] and [5] provide a 
few examples of PDFs derived for stars in different positions in 
the CMD. They illustrate the typical cases of PDFs with single, 
well-defined peaks (Fig.|4} as well as some cases for which the 
parameters cannot uniquely be determined (Fig.|5|i. 

Two cases for which the parameter determination works ex- 
cellently are HD 62644 and HD 34642 (Fig.g] upper half). The 
PDFs present a single prominent peak and a modest dispersion 
for all estimated parameters, permitting a clear identification 
of the most likely parameter value by the mean and its uncer- 
tainty as defined by the standard error. The errors in determin- 
ing mass and age are however much smaller for HD 62644 than 
for HD 34642, even if these two stars have a similar (and very 
small) parallax error. This is so because HD 62644 is still a sub- 
giant, whereas HD 34642 is an authentic red giant. The evolu- 
tionary tracks for different masses for HD 62644 run horizon- 
tally and well separated in the CMD, whereas for HD 34642 
they are more vertical and much closer. The same input error 
bars will in general produce larger output errors for a giant than 
for a subgiant. 

One notices that the errors in parameter determination 
for these two cases originate primarily from the error in the 
Hipparcos parallax, and are just slightly increased by the ef- 
fect of a 0.05 dex error in [Fe/H]. In fact, the [Fe/H] error 
involves the age-metallicity degeneracy, especially for the red 
giant HD 34642. It is also worth noticing that parameters like 
g, r e ff, and R depend slightly on the measured metallicity and 
its error. Acceptable values for these parameters could have 



been derived assuming a fixed - say equal to solar - metal- 
licity. However, our method has the advantage of accounting 
fully for the subtle variations of the colour and bolometric cor- 
rection scales with metallicity, thereby slightly reducing the fi- 
nal errors in the parameter determination when the observed 
[Fe/H] is taken into account. In contrast, for estimating t and 
M, a proper evaluation of the metallicity and its error turns out 
to be absolutely necessary. 

Other well-behaved cases are HD 50778 and HD 125560 
(Fig. EJ, which are located in the upper part of the RGB, one 
well above and the other close to the red clump region. A robust 
mass and age determination is possible for these stars, but their 
mass PDFs show a faint secondary peak close to a prominent 
primary peak. For HD 50778, the primary peak corresponds to 
~ 0.95 M for a star in the phase of first-ascent RGB, whereas 
the secondary peak corresponds to the early-AGB phase of a 
0.7 M star. There is no way to distinguish a priori between 
these two evolutionary phases from our observations, but for- 
tunately the first-ascent RGB case turns out to be much more 
likely (due to its longer evolutionary timescale) than the early- 
AGB phase. In the case of HD 125560, the primary peak in 
the PDF corresponds to a red clump (core He-burning) phase 
of a 1.1 M star, whereas the secondary peak corresponds to a 
1.6 M star in the first-ascent RGB phase. Similar results, with 
the presence of small secondary peaks in the PDF, are common 
for stars in the upper part of the CMD. 

Parameter estimation is much more difficult for stars like 
HD 1 1977 and HD 174295 (Fig.|5). These stars are in the mid- 
dle of the most degenerate region of the CMD, namely in the 
"loop" region of red clump stars of different masses. As illus- 
trated in Girardi et al. (1998, their figure 1), core He-burning 
stars with the same metallicity and different mass form a com- 
pact loop which starts for low masses in the blue end as an ex- 
tension of the horizontal branch, reaches its reddest colour as 
the mass increases and then turns back into the blue direction. 
The luminosity along the same mass sequence first increases 
slowly by some tenths of a magnitude, decreases sharply by 
0.5 mag at about 2 M , and turns towards much higher lumi- 
nosities as the mass increases further. Such a complex pattern 
in the CMD implies that for some stars in the middle of such 
loops, two mass (and age) values may become similarly likely, 
resulting in bi-modal PDFs. Moreover, the shapes of such PDFs 
may become sensitive to even small changes in [Fe/H], thus 
increasing further the uncertainty of the determination of the 
stellar parameters. Ill-behaved cases like the ones illustrated in 
Fig.|5]account for less than a fifth of our sample. 

3.3. Results and checks 

Results for all our sample stars are presented in columns 6 to 
17 of Table |6] Except for (B-V)o, we have determined mean 
values and standard errors using PDFs of logarithmic quanti- 
ties, and subsequently converted these values into linear scales. 
We have done so because linear quantities (age, mass, radii, 
etc.) are more commonly used and are considered to be more 
intuitive than the corresponding logarithms. Whenever possi- 
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Fig. 4. Examples of probability distribution functions (PDFs), illustrating the majority of well-behaved cases in our sample for 
which good mass and age estimates are possible. For each star, one panel presents the position in the HR diagram (dot). The five 
remaining panels show the PDF for log f, M, \ogg, log/?, and (B-V). The dotted lines show the PDF assuming no error in the 
[Fe/H] determination; in this case, the PDF width is mostly determined by the parallax errors. The solid lines show the slightly 
broader PDF obtained by assuming an internal cr [Fe / H ] of 0.05 dex. HD 34642 and HD 62644 (upper panels) represent stars in 
the lower part of the RGB for which single and well defined peaks in the PDF are typical. HD 50778 and HD 125560 (bottom 
panels), being located in the upper part of the RGB, exhibit small secondary peaks in the mass and age PDFs, which represent a 
small probability that the star belongs to an evolutionary phase different from the one causing the main peak. 



ble, however, we will use the original error bars obtained with 
the logarithmic scale. 



We have applied a few checks to what extent our method 
for parameter estimation is reliable. 



3.3.1 . "Colour excess" (B-V) - (B-V) 

The first test involves a comparison of the estimated intrinsic 
colours (B-V)o with the observed colours (B-V). Their differ- 
ences are presented in Figs. |6] and as a function of distance 
and r e ff. Since the (B-V)o values result essentially from the 
spectroscopic T^, error bars in the individual (B-V)-(B-V)o 
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Fig. 5. Same as Fig. 0] illustrating two cases for which the estimate of stellar parameters is not unique. Both stars are in the red 
clump region of the HR diagram. Their parameters can be well reproduced by clump stars of different masses/ages, which are 
represented by the two main peaks in their PDFs. 




-0.2 



100 150 
distance (pc) 




-0.2 



Fig. 6. Difference between the observed (B—V) values, and 
the (B-V\) values derived from the PDF method, as a function 
of distance from the Sun. Notice the absence of a clear trend 
of (B—V) - (B-V)o with distance. The data points shown as 
crosses denote the cases which we consider as "outliers" and 
which are excluded from our statistical analysis. 



values reflect mostly the 70 K error assumed for T e g, except 
for the few stars for which the (B—V) error of the Hipparcos 
catalogue was significant (see also Hig . |3]> . 



Fig. 7. The same as Fig. [6] but as a function of effective tem- 
perature. Notice that the small differences found for most of the 
stars (less than 0.05 mag) appears to be a function of T e s. 



As can be clearly seen in Fig. [6] there is no marked in- 
crease of (B-V) - (B-V)o with distance, as can be expected 
for a sample with distances less than ~ 200 pc and correspond- 
ingly little reddening. Taking the diffuse interstellar absorption 
of dAy/dr = 0.75 mag/kpc (Lynga 1982), one expects a color 
excess of just E&-v = 0.05 mag at 200 pc. This expectation is 
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consistent with our (B-V)-(B — V) data, although comparable 
with their dispersion. 

Another aspect to notice in Fig.[6]is the small dispersion of 
the (B - V)—(B - V)o values we obtained for the bulk of the sam- 
ple stars. If we disregard two outliers with \(B -V) — (B— V)o\ ^ 
0.1, we find an unweighted mean of (B-V)-(B-V) = -0.009 
with a scatter of 0.03 1 . This scatter (excluding outliers) can be 
considered the typical error of our PDF method for determining 
the intrinsic colour of our sample. 

Fig.Qshows how (B— V)— (B— V)q depends on T e g. There is 
evidently a correlation between (B—V)—(B— V)o and T e ff , with 
a minimum difference at ~ 4700 K and a maximum difference 
at ~ 4000 K. It is likely that this correlation is caused by errors 
in the theoretical r e jf-colour relation adopted in the Girardi et 
al. (2002) isochrones, which amount to less than 0.05 mag, or 
equivalently to about 100 K for a given (B-V)o l . In addition, 
even if our T e s -colour scale were perfectly good, (B—V) starts 
to become intrinsically a poor T e ff indicator for the coolest gi- 
ants. 

Notice that the possible systematic errors of 0.05 mag in 
our adopted r e ff-colour relation would imply errors smaller 
than 0.03 for the V-band bolometric corrections adopted for the 
same isochrones, which would then be the maximum mismatch 
between theoretical and observational My values. We conclude 
that these errors are small enough to be neglected in the present 
work. 

The two outliers with high \(B-V) - (B-V\)\ in Fig. |6] 
(HD 22663 and HD 99167) can be explained as (1) stars with a 
significant reddening; (2) stars for which Hipparcos catalogue 
has a wrong entry for (B - V); or (3) stars for which our parame- 
ter estimation (including r e ff, [Fe/H], and/or (B-V)o) substan- 
tially failed. We consider the third alternative as being the most 
likely one, and hence exclude these two stars from any of the 
statistical considerations that follow in this paper. 

3.3.2. Surface gravities 

Another important check is the comparison between our esti- 
mated log g values with those derived independently from spec- 
troscopy (see Sect. |3, presented in Fig.|S] The PDF-estimated 
values tend to be systematically lower than the spectroscopi- 
cally derived ones. Again ignoring the two outliers of the pre- 
vious section, the mean difference is -0.20 dex with a stan- 
dard deviation of 0.14 dex. Such an offset would indicate, for 
instance, that our method underestimates stellar masses by a 
factor of about 1 .6, which however can be excluded given our 
results for the two Hyades giants (see below). An alternative 
explanation is that the spectroscopic \ogg values are simply 
too high. 

The latter interpretation is supported by the consideration 
that gravity is determined by imposing ionization balance; this 
means that the abundance found for the nine Fe II lines is the 
same as the one retrieved for the (more than 70) Fe I lines. 
This procedure implies that spectroscopic gravities depend, in 
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A small contribution to this behaviour with r eff may be due to 
modest reddening of the most distant stars, which have a tendency 
towards brighter absolute magnitudes and lower Tgg, 



2 3 
spectroscopic log(g) 



Fig. 8. Comparison between the log g values derived from the 
spectra (abcissa) and from the photometry by means of our 
PDF method (ordinate). The solid line represents a one-to-one 
relation. The upper panel shows the differences between esti- 
mated and spectrocopic values. 



addition to the adopted line oscillator strengths, to the interplay 
between the stellar parameters in the derivation of abundances. 
This can be fairly complex in Pop I giants, where the Fe I vs. 
Fe II abundance depends not only on gravity, but also quite 
strongly on effective temperature. As an example, see Pasquini 
et al. (2004, their table 4), where the dependence of Fe I and 
Fe II on log g, T e g and £ is analyzed for one Pop I giant and 
the same set of lines. A systematic shift of 100 K in T e g would, 
for instance, produce a 0.2 dex shift in logg without changing 
substantially the derived Fe abundance. 

The disagreement between the gravity values obtained from 
spectroscopy and from parallaxes has been known for a long 
time (e.g., da Silva 1986), and the problem did not disap- 
pear despite the improvements of models and parallax mea- 
surements. Nilsen et al. (1997) compare the Hipparcos-based 
gravities with the values obtained from spectroscopy by sev- 
eral authors and conclude that differences between the two 
methods could become larger than a factor two (0.3 in \ogg 
). This can have various causes, like non-LTE effects on Fe I 
abundances, or thermal inhomogeneities. We conclude that our 
spectroscopic gravities, like that of other authors, are systemat- 
ically overestimated. Note also that Monaco et al. (2005, their 
figure 6) find that spectroscopic log g values correlate with mi- 
croturbulent velocities in such a way that a systematic er- 
ror of just 0.07 km s _1 in £ would be sufficient to cause the 
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Fig. 9. Comparison of apparent stellar diameters 9 from the 
CHARM2 catalogue (abcissa) and derived photometrically 
with our PDF method (ordinate). The entries in CHARM2 were 
corrected by limb darkening as detailed in the text. The solid 
line represents the one-to-one relation. The large dots with 
HD identification indicate stars with 9 measurements obtained 
via lunar occultation and/or LBI techniques. The upper panel 
shows the relative differences A9/theta, i.e., the differences be- 
tween estimated and measured values, divided by the measured 
6. 



-0.20 dex offset which we find in log g. We point out that the 
-0.20 dex offset in \ogg would have a negligible effect on our 
derived [Fe/H] values, that are mostly based on the gravity- 
insensitive Fe I lines (see section 3). 

3.3.3. Stellar radii and apparent stellar diameters 

Another check regards the stellar radii R, which can be eas- 
ily converted into apparent stellar diameters 6 using Hipparcos 
parallaxes, and then compared to observations. The observed 
apparent diameters 9 are taken from the CHARM2 catalogue 
(Richichi et al. 2005), which in most cases consist of indi- 
rect estimates of stellar diameters using fits to spectropho- 
tometric data, and are available for about one third of our 
sample. "Direct" diameter determinations are available only 
for five of our giants, based on lunar occultation (HD 27371 
and HD 27697; the Hyades giants) and on long baseline in- 
terferometry (LBI; HD 27697, HD 81797, HD 113226, and 
HD 189319). In cases for which more than one 9 determination 
was available, we used either the most accurate one (i.e. the 



one with a substantially smaller error) or the most recent one 
when tabelled errors were similar. Finally, we have corrected 
all uniform-disk measurements by limb darkening (LD) using 
the extensive tables provided by Davis et al. (2000), that are 
based on Kurucz (1993ab) model atmospheres. For the 20 sam- 
ple stars with LD-corrected diameters in CHARM2, our correc- 
tions agree perfectly with those provided there. The mean LD- 
correction for these giants is 3.7 + 3.3 %, which is well below 
the ~ 12 % lcr relative error of our individual 9 estimates. 

The comparison between our derived 9 values and the 
LD-corrected CHARM2 diameters is presented in Fig. |9| 
Here again, the comparison is very satisfactory, with the es- 
timated minus observed difference -0.21 mas with a scatter of 
0.32 mas. Alternatively, we looked at the fractional differences 
A#/#charm2 (upper panel of Fig.|5}: its unweighted mean value 
is of 69/9 = -0.06 with a r.m.s. scatter of 0.06. This scatter also 
is well below the fractional error of individual 9 estimates, of 
erg/9 =i 0.12 (see the error bars in the upper panel of Fig. [9}, 
which are largely due to Hipparcos parallax errors. 

A similar level of agreement is obtained for stars with "di- 
rect" 9 mesurements. For all the other stars in Fig. |5] the "in- 
direct estimates" based on fits to spectrophotometric data (and 
especially infrared data) presented in CHARM2 correlate very 
well with our values. This is a remarkable result. It indicates 
that by using just two visual passbands (B and V as in the 
present work) for giants of known distance, T e ff and metallicity, 
it is possible to obtain diameter estimates of a quality similar 
to that obtained by more sophisticated methods based on multi- 
band spectrophotometry. 

Our results indicate a very successful and robust estimation 
of the stellar parameters (B - V)o, log g, and R, for the bulk of 
our sample stars. 

3.3.4. Ages and masses 

With respect to the parameters t and M, the few checks at our 
disposal address the correlations with [Fe/H], and the results 
for the two Hyades giants in our sample. Fig. [^presents the 
mass-metalicity plot. It shows a clear pattern: low-metallicty 
giants (with [Fe/H] < -0.4) are present only among the stars 
with the lowest masses. All stars with M > 1.2 M e are char- 
acterised by a mean solar metallicity (0.00 dex) and a small 
standard deviation of 0. 12 dex. The same data, when presented 
as an age-metallicity plot in Fig.[^ indicates a more scattered 
pattern but with the similar indication of metal-poor stars being 
present with an age above ~ 5 Gyr. The present data points to 
an age-metallicity relation with a large scatter for the highest 
ages. This scatter is consistent with other results in the litera- 
ture, at least for the highest ages (see Sect.|4}. 

Our sample contains two Hyades giants for which good- 
quality age and mass estimates are available: HD 27371 with 
t = 0.53 ± 0.09 Gyr, M = 2.70 + 0.13 M , and HD 27697 
with t = 0.67 ±0.13 Gyr, M = 2.54 + 0.14 M . The Hyades 
turn-off age, as derived from models with overshooting, is 
0.625 + 0.05 Gyr (Perryman et al. 1998). 2 This value is con- 

2 Ages close to 0.63 Gyr are also indicated by the Hyades binaries 
V 818 Tauri, 51 Tauri, and 8 2 Tauri (Lastennet et al. 1999, in their 
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Fig. 10. Mass-metallicity relation for our sample stars. 
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Fig. 11. Age-metallicity relation for our sample stars. 



sistent to within lcr with our estimated ages for HD 27697 and 
HD 27371. 

Although the error in our age estimates for individual stars 
is unconfortably large compared to the typical error of clus- 
ter turn-off ages, the pair of Hyades giants provides the best 
evidence that the J0rgensen & Lindegren (2005) method works 

table 6), whereas the white dwarfs provide a lower age limit of 0.3 ± 
0.03 Gyr (Weidemann et al. 1992). 



well for estimating the ages of giants. Unfortunately, additional 
checks using other clusters are apparently impossible at this 
moment. Although other well-studied clusters with excellent 
turn-off ages exist (e.g. M67 and Praesepe), they do not belong 
to our sample and have significantly smaller parallaxes. 

To summarize, we conclude that our method has provided 
excellent determinations of stellar parameters, especially for 
(B-V)o, logg, R and for which the uncertainties of the 
method were intrinsically low (with a few exceptions), and has 
been confirmed by independent data. For the stellar ages and 
masses, however, our determinations turn out to be intrinsically 
more uncertain, as demonstrated by the larger error bars we ob- 
tained. Although we have some indication from the Hyades gi- 
ants that our age scale is not very inaccurate, another indepen- 
dent check with other mass and age data would be desireable. 

4. More about the age-metallicity relation 

The Solar Neighbourhood age-metallicity relation (AMR) pro- 
vides basic information about the chemical evolution of the 
Milky Way's disk with time, and has for long been used to con- 
strain evolutionary models of our Galaxy. We refer to Carraro 
et al. (1998), Feltzing et al. (2001) and Nordstrom et al. (2004) 
for recent determinations of the AMR and a general discussion 
of its properties. 

Since we have derived the AMR from a completely new 
sample and use a relatively new method, it is important to il- 
lustrate how our results compare with other determinations. 
Note also that our sample is appropriate to do this because the 
stars were not chosen using any criterion of age, abundance or 
galatic velocity. 

First, we transform our AMR data of Fig.^^to a more sim- 
ple function of age. To do so, for each age bin Af, we determine 
its cumulative [Fe/H] by adding the measured [Fe/H] of each 
star weighted by its probability of belonging to Af . This prob- 
ability is given by the age PDF of each star, integrated over 
the Af interval. We obtain a cumulative [Fe/H] distribution for 
each age bin, from which we derive the mean [Fe/H] value and 
dispersion. Since the metallicity value of a single star is spread 
over several age bins (just as for the age PDF), the effective 
number of stars per bin, (AO, may be a fractional number. We 
choose age bins wide enough to provide (N) > 1 for all ages. 
Note that the data points obtained this way for different age 
bins are not independent. The results are presented in Fig. El 
and Table [5] where the mean metallicity and its dispersion is 
shown as a function of age. 

Our results present a few notable characteristics: 

1. The AMR relatively flat up to the largest ages. From its 
present-day solar's value, the mean metallicity has dropped 
to just -0.2 dex at an age of 12 Gyr. This result appears to 
be in qualitative agreement with Nordstrom et al. (2004). 
Other authors arrive at a somewhat lower mean metallicity 
at the oldest ages (e.g., Carraro et al., 1998; Rocha-Pinto et 
al., 2000; Reddy et al., 2003). As discussed by Nordstrom 
et al. (2004) and Pont & Eyer (2004), selection against old 
metal-rich dwarfs may have contributed to the steeper de- 
cline of [Fe/H] with age found in many of the previous stud- 
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Fig. 12. Age-metallicity relation for 1 Gyr age bins (see 
Tabled. 

Table 5. Age-metallicity relation derived from our sample for 
1-Gyr wide age bins. 



{t) (Gyr) 


<[Fe/H]> 


C [Fe/H] 


(N) 


0.5 


0.01 


0.09 


17.1 


1.5 


-0.04 


0.15 


10.5 


2.5 


-0.03 


0.17 


6.1 


3.5 


0.00 


0.18 


9.1 


4.5 


-0.02 


0.20 


5.3 


5.5 


-0.04 


0.22 


5.3 


6.5 


-0.08 


0.24 


4.7 


7.5 


-0.10 


0.26 


3.4 


8.5 


-0.12 


0.27 


2.3 


9.5 


-0.15 


0.29 


2.5 


10.5 


-0.23 


0.30 


4.2 


11.5 


-0.28 


0.30 


1.6 



ies. It is likely that such selection effects are absent in our 
sample of giants. 

2. The metallicity dispersion tends to increase with age, so 
that one finds a large [Fe/H] dispersion of about 0.3 dex 
(lcr) for the highest ages. Part of this trend may result from 
the increase in the age error with age (evident in Fig.llH. 
which causes stars of very different ages to contribute to the 
mean [Fe/H] of the same age bins. It seems clear from our 
data that for ages larger than about 4 Gyr, the [Fe/H] dis- 
persion becomes considerably larger than the observational 
errors. This is in fairly good agreement with the results of 
other authors. 

3. For all ages lower than about 4 Gyr, we find that the [Fe/H] 
dispersion is comparable to the typical observational error. 
The most striking result hoewever is given by the to 1- 
Gyr age bin, which is very well populated (17.1 stars) and 
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typically very low so that the confu- 
sion witn otner age bins is practically absent. In this case, 
the [Fe/H] dispersion is 0.09 dex, and compares well with 
the observational [Fe/H] errors of 0.05 dex (internal) and 
0.1 dex (external). This result is in contradiction into the 
claims by most authors (Carraro et al., 1998; Feltzing et al., 
2001; Nordstrom et al., 2004, etc), who find large [Fe/H] 
spread at all ages. 

Why do our results differ from other authors, in particular 
for the youngest stars? The main difference is likely to result 
from the different kinds of stars investigated in the mentioned 
above studies. We study only giants whereas most authors use 
field dwarfs and a few subgiants. As we have demonstrated, 
the age determination for the youngest giants appears to be 
quite reliable. The same may not be true for dwarfs which, be- 
ing located close to the main sequence, are in a more degener- 
ate region of the HR diagram. The results of Edvardsson et al 
(1993) and Nordstrom et al. (2004), which represent the best 
age estimates for a limited section of the main sequence, report 
considerable uncertainities for the ages of most of their stars. 
In addition, their [Fe/H] determination is based on ubvy/3 pho- 
tometry, which further limits their sample to an even smaller 
section of the main sequence. These problems do not exist 
for giants with spectroscopic [Fe/H] and T e fi determinations. 
Reliable measurements can be obtained all along the RGB, and 
the worst problems appear not to be related to selection effects, 
but to the intrinsic age errors illustrated in Sect. 13.21 We con- 
clude that giants may be better targets for the study of the Solar 
Neighbourhood AMR than dwarfs. As an alternative, subgiants 
may be even better because they always provide reliable age 
estimates, as illustrated in Fig. |4] and as already explored by 
Thoren et al. (2004). 

The interpretation of our AMR result indicates that stars 
in the Solar neighbourhood are formed from interstellar matter 
of quite homogeneous chemical composition. As we observe 
older stars, we start sampling stars born in different Galactic 
locations 3 , and hence we see a more complex mixture of chem- 
ical composition. 

5. Discussion 

5.1. RV variation 

Setiawan et al. (2004) have studied the trend of radial velocity 
(RV) variability - determined by the standard deviation cr RV - 
with My, detecting an increase of variability along the RGB. 
We can now do the same for other stellar parameters. FigurelT3l 
shows our star sample in the HR diagram obtained from the 
information of Tabled We note that the cr RV represents the to- 
tal standard deviation without regard to the timescales involved 
(i.e. short- versus long term variations). The 1 1 binaries of the 
sample which are reported by Setiawan et al. (2004) are ex- 
cluded from this plot. There is a trend of increasing cr RV with 
position along the RGB, which becomes more evident if we 

3 Different Galactic locations mean different radial positions in the 
thin disk. The oldest age bins may also contain a few thick disk and 
halo stars. 
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Fig. 13. HR diagram of our star sample, excluding 1 1 bina- 
ries. Stellar parameters were taken from Table [6] Each cir- 
cle diameter is proportional to cr RV , ranging from ~ 30 to 
280 ms 4 . Stars that host substellar companions are marked 
with full squares for suspected brown dwarfs and full triangles 
for suspected giant planets. 

exclude the stars that host sub-stellar companion candidates 
(marked with different symbols). 

A similar trend is seen in the plot of <x RV against stellar 
radius in Fig.^J The two stars suspected to host brown dwarf 
companions, HD 27256 and HD 224533, stand out as showing 
too large a cr RV for their radii. If these two stars are excluded, 
an unweighted least squares fit to the data results in the relation 

o-RvCms- 1 ) = (1.76 ±03l)R(R Q ) + (38.1 ±6.0). 

The Spearman rank correlation coefficient is 0.45. This result, 
although supporting the increase of cr RV with R, is not signifi- 
cant because <x RV is limited by the long-term precision of our 
RV measurements (about 25 ms 4 , see Setiawan et al. 2003, 
2004). 

5.2. Distribution of metallicity 

Another major result of this study is the derivation of the age- 
metallicity relation for the Solar Neighbourhood (see Fie. II II . 
Its behaviour agrees, in general, with most previous determi- 
nations in the literature, except for the very low [Fe/H] spread 
that we find for the youngest ages, which is comparable to the 
observational error. The main novelty with respect to previous 
results is that we derive the AMR using data for field red gi- 
ants only, whereas the majority of present-day determinations 
have used samples of field dwarfs (including just a small frac- 
tion of subgiants), or giants belonging to open clusters. It worth 
noticing that we used very simple data in our determinations, 



Fig. 14. <t rv versus stellar radius diagram for our sample, ex- 
cluding 11 binaries. Refer to Fig. [D] for the meaning of the 
symbols. 

namely the BV photometry together with Hipparcos parallaxes 
and measured values for T e ff and [Fe/H] . 

Of course, the same work can be extended to all giants in 
Hipparcos catalog, once we have obtained homogeneous T e g 
and [Fe/H] determinations for them. This opens the possibility 
of improving considerably the statistics and reliability of the lo- 
cal age-metallicity relation, simply by acquiring spectroscopic 
data for a larger sample of bright giants, and performing the 
same abundance and parameter analysis as in the present work. 

Moreover, similar methods can be applied to nearby galax- 
ies with well-known distances, once we have available both the 
photometry and spectroscopy for a sufficiently large number of 
their red giants. Zaggia et al. (in preparation) use a procedure 
similar to ours to derive the age-metallicity relation of the Sgr 
dSph galaxy. 

5.3. Stars hosting low mass companions and planets 

In the course of earlier studies we have identified three stars 
as candidates to host low mass companions: two with planets 
(HD 47536 and HD 122430) and one with either a brown dwarf 
or a planet companion (HD 1 1977). We investigate a range of 
stellar masses larger than the range of masses usually investi- 
gated with radial velocity techniques. Not only can we provide 
better constraints on the companion mass, but we can also in- 
vestigate to what extent the conditions for companion forma- 
tion differ within the mass range. Only one of our three stars, 
HD 122430, has nearly solar metallicity ([Fe/H] = -0.05), 
while HD 11977 is slightly sub-solar ([Fe/H] = -0.21) and 
HD 47536 ([Fe/H] = -0.68) is the most metal poor star of the 
sample. Schuler et al. (2005) derived a metallicity of [Fe/H] = 
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-0.58 for the K giant star hosting planet HD 13189. Although 
we are considering a small number of objects, this result seems 
at odds with what has been found for dwarf stars hosting gi- 
ant exoplanets, which are preferentially metal rich (e.g. Santos 

2004) . However, most RV planet search programs have concen- 
trated on solar mass stars and two of our planet hosting giant 
stars have masses considerably larger than solar. In the case of 
HD 13189 the host star has a mass of 3.5 M Q (Schuler et al. 

2005) . At the present time is is unknown what role stellar mass 
plays in the process of planet formation and for massive stars 
this may be a more dominant factor than the metallicity. Any in- 
vestigation of the metallicity-planet relation among giant stars 
should focus on those in the same mass range. It may be that 
for a given mass range stars with higher metal abundances still 
tend to host a higher frequency of giant planets. 
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Table 6. Stellar parameters as derived from the observed My, T tf! , and [Fe/H] via the PDF method. Errors in T eS and [Fe/H] are 70 K and 
0.05 dex for all stars. The [Fe/H] values shown here were corrected by a zero-point offset of -0.07 dex. E^ v stands for the observed (B- V) 
minus the estimated intrinsic (B-V)o- 



HD 




M v 
(mag) 


T e ff 
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[Fe/H] 
(dex) 




t 
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M 

(M ) 


(c.g.s.) 


R 

(Re) 
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e 
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2114 


-0.53 


+ 


0.41 


5288 


-0.03 


0.32 


+ 


0.07 


3.16 


± 


0.29 


2.62 + 0.18 


13.8 + 3.5 


-0.00 
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0.04 


0.71 


± 


0.31 


2151 


3.45 
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0.01 
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-0.03 


5.12 
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1.03 


1.17 
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0.05 


4.02 + 0.04 
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0.01 
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0.07 
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0.95 
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0.12 
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4084 


-0.29 


8.11 


+ 


2.37 


1.03 


+ 


0.10 


1.46 + 0.10 


30.1 ±2.9 


0.03 


± 


0.05 


3.62 


± 


0.60 


61935 


0.71 


± 


0.08 


4879 


-0.01 


1.18 


+ 


0.42 


2.02 


+ 


0.29 


2.71+0.09 


10.1 ±0.5 


0.00 


± 


0.04 


2.11 


± 


0.18 


62644 


3.13 


± 


0.04 


5526 


0.12 


3.78 


+ 


0.19 


1.31 


± 


0.03 


3.78 + 0.03 


2.35 ± 0.09 


-0.00 


± 


0.03 


0.90 


± 


0.05 


62902 


1.17 


+ 


0.15 


4311 


0.33 


6.95 


+ 


2.57 


1.10 


+ 


0.16 


2.29 ± 0.09 


11.9 ±1.1 


0.06 


± 


0.04 


1.52 


± 


0.24 


63697 


0.72 


± 


0.12 


4322 


0.13 


4.62 


± 


2.17 


1.26 


± 


0.21 


2.16 + 0.10 


14.8 ± 1.3 


-0.02 


± 


0.04 


1.77 


± 


0.26 


65695 


0.51 


± 


0.16 


4468 


-0.14 


4.29 


± 


2.17 


1.27 


± 


0.20 


2.19 + 0.10 


14.4 ±1.4 


0.01 


± 


0.04 


1.75 


± 


0.30 


70982 


0.34 


± 


0.15 


5089 


-0.03 


0.61 


+ 


0.12 


2.54 


± 


0.15 


2.72 ± 0.05 


11.0 ± 0.8 


-0.03 


± 


0.02 


0.72 


± 


0.10 


72650 


0.45 


+ 


0.15 


4310 


0.06 


4.05 


+ 


1.94 


1.32 


+ 


0.20 


2.06 + 0.11 


17.0 ± 1.6 


-0.01 


± 


0.04 


1.05 


± 


0.18 


76376 


0.17 


+ 


0.20 


4282 


-0.10 


5.39 


+ 


2.56 


1.20 


± 


0.17 


1.87 + 0.13 


20.1 ±2.4 


0.02 


± 


0.04 


1.44 


± 


0.30 


81797 


-1.69 


+ 


0.09 


4186 


0.00 


0.42 


+ 


0.16 


3.03 


± 


0.36 


1.48 + 0.10 


50.5 ±4.0 


0.04 


± 


0.04 


8.65 


± 


1.06 


83441 


0.89 


+ 


0.15 


4649 


0.10 


2.98 


+ 


1.57 


1.40 


+ 


0.28 


2.48 + 0.11 


10.8 ±0.6 


-0.03 


± 


0.03 


0.98 


± 


0.12 


85035 


2.62 


+ 


0.14 


4680 


0.12 


6.42 


± 


2.22 


1.17 


+ 


0.11 


3.12 + 0.08 


4.77 ± 0.37 


-0.03 


± 


0.03 


0.58 


± 


0.08 


90957 


0.21 


± 


0.17 


4172 


0.05 


5.57 


± 


2.46 


1.18 


+ 


0.18 


1.79 + 0.13 


22.1 ±2.3 


0.02 


± 


0.04 


1.73 


± 


0.32 


92588 


3.57 


± 


0.06 


5136 


0.07 


4.72 


± 


0.53 


1.24 


± 


0.04 


3.75 + 0.03 


2.36 ±0.12 


-0.03 


± 


0.03 


0.64 


± 


0.05 


93257 


1.81 


± 


0.09 


4607 


0.13 


4.95 


± 


2.09 


1.26 


± 


0.15 


2.75 ± 0.09 


7.51 ±0.51 


-0.03 


± 


0.03 


1.28 


± 


0.14 


93773 


0.72 


± 


0.21 


4985 


-0.07 


0.96 


± 


0.27 


2.15 


± 


0.21 


2.79 ± 0.07 


9.44 ± 0.76 


-0.03 


± 


0.02 


0.61 


± 


0.11 


99167 


-0.43 


+ 


0.18 


4010 


-0.36 


8.92 


± 


2.10 


0.99 


± 


0.08 


1.33 + 0.10 


34.1 ±3.8 


0.14 


± 


0.04 


2.84 


± 


0.55 


101321 


1.72 


± 


0.19 


4803 


-0.14 


3.96 


± 


1.62 


1.31 


± 


0.14 


2.88 + 0.10 


6.62 ± 0.66 


-0.05 


± 


0.03 


0.59 


± 


0.11 


107446 


-0.63 


± 


0.09 


4148 


-0.10 


2.83 


+ 


1.20 


1.42 


± 


0.22 


1.52 + 0.11 


32.9 ±2.5 


-0.02 


± 


0.07 


4.38 


± 


0.51 


110014 


-0.29 


+ 


0.19 


4445 


0.19 


0.86 


+ 


0.34 


2.30 


+ 


0.28 


2.06 + 0.10 


22.6 ± 2.2 


-0.05 


± 


0.04 


2.15 


± 


0.39 


111884 


0.58 


+ 


0.17 


4271 


-0.06 


7.27 


± 


2.62 


1.12 


+ 


0.11 


2.02 + 0.10 


16.6 ±1.7 


0.01 


+ 


0.04 


1.33 


± 


0.24 


113226 


0.37 


+ 


0.06 


5086 


0.09 


0.56 


+ 


0.09 


2.64 


+ 


0.11 


2.71+0.04 


11.4 ±0.5 


-0.07 


± 


0.03 


3.38 


± 


0.23 


115478 


0.53 


± 


0.20 


4250 


0.03 


6.05 


+ 


2.54 


1.18 


+ 


0.14 


2.00 + 0.11 


17.4 ±2.0 


-0.03 


± 


0.04 


1.77 


± 


0.36 


122430 


-0.15 


+ 


0.21 


4300 


-0.05 


3.11 


± 


1.70 


1.39 


± 


0.27 


1.83 + 0.12 


22.9 ±2.3 


0.01 


± 


0.04 


1.60 


± 


0.32 


124882 


-0.35 


± 


0.09 


4293 


-0.24 


4.30 


+ 


2.07 


1.21 


+ 


0.21 


1.69 + 0.11 


24.9 ± 1.8 


0.01 


± 


0.30 


2.71 


± 


0.31 


125560 


1.01 


± 


0.09 


4472 


0.16 


5.21 


+ 


2.28 


1.14 


± 


0.19 


2.36 ± 0.08 


11.3 ±0.5 


0.01 


± 


0.03 


1.79 


± 


0.15 


131109 


-0.22 


± 


0.16 


4158 


-0.07 


5.14 


± 


2.43 


1.18 


± 


0.20 


1.64 + 0.11 


26.3 ±2.6 


0.05 


± 


0.04 


1.87 


± 


0.32 


136014 


0.83 


+ 


0.21 


4869 


-0.46 


5.32 


+ 


2.80 


1.05 


± 


0.21 


2.44 ± 0.09 


9.86 ± 0.66 


0.01 


± 


0.03 


0.78 


± 


0.13 
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HD 




My 

(mag) 


(K) 


[Fe/H] 

(dex) 




t 

(Gyr) 




M 

(M ) 




(c.g.s.) 


R 

(Re) 




Eb-v 
(mag) 




6 

(mas) 


148760 


1.87 


+ 


0.12 


4654 


0.13 


4.16 


+ 


1.61 


1.32 


+ 


0.15 


2.84 


± 


0.09 


6.99 + 0.54 


-0.05 


± 


0.03 


0.94 


+ 


0.13 


151249 


-1.14 


± 


0.14 


3886 


-0.37 


7.11 


+ 


2.34 


1.02 


+ 


0.12 


0.92 


± 


0.12 


55.9 + 7.3 


0.02 


+ 


0.04 


5.42 


± 


1.06 


152334 


0.30 


+ 


0.09 


4169 


0.06 


5.80 


+ 


2.26 


1.19 


+ 


0.14 


1.84 


± 


0.10 


21.0+1.6 


0.00 


± 


0.04 


4.23 


+ 


0.49 


152980 


-0.79 


+ 


0.16 


4176 


0.01 


1.70 


± 


0.57 


1.74 


± 


0.24 


1.59 


± 


0.11 


33.7 + 3.4 


0.04 


+ 


0.04 


3.36 


+ 


0.58 


159194 


1.61 


+ 


0.20 


4444 


0.14 


6.76 


+ 


2.42 


1.15 


+ 


0.13 


2.56 


± 


0.11 


8.97 + 0.98 


0.02 


+ 


0.03 


0.78 


+ 


0.16 


165760 


0.33 


± 


0.13 


5005 


0.02 


0.60 


+ 


0.12 


2.59 


± 


0.15 


2.70 


± 


0.05 


11.5+0.7 


-0.04 


± 


0.04 


1.46 


+ 


0.18 


169370 


1.61 


± 


0.18 


4460 


-0.17 


9.31 


± 


1.90 


1.05 


± 


0.06 


2.56 


± 


0.07 


8.61+0.70 


0.03 


± 


0.02 


0.93 


± 


0.15 


174295 


0.41 


± 


0.14 


4893 


-0.24 


1.87 


± 


0.86 


1.60 


± 


0.31 


2.54 


± 


0.10 


10.8+0.7 


-0.02 


± 


0.03 


1.12 


± 


0.14 


175751 


0.82 


± 


0.12 


4710 


0.01 


2.78 


+ 


1.39 


1.42 


+ 


0.27 


2.50 


± 


0.10 


10.7 + 0.5 


-0.05 


± 


0.06 


1.56 


± 


0.16 


177389 


2.48 


+ 


0.05 


5131 


0.02 


1.87 


± 


0.15 


1.65 


+ 


0.05 


3.43 


± 


0.04 


3.93 + 0.19 


-0.01 


± 


0.03 


0.99 


± 


0.07 


179799 


2.35 


+ 


0.17 


4865 


0.03 


3.61 


+ 


1.21 


1.36 


+ 


0.13 


3.17 


± 


0.09 


4.80 ± 0.43 


-0.05 


± 


0.03 


0.65 


+ 


0.11 


187195 


1.29 


+ 


0.14 


4444 


0.13 


6.44 


+ 


2.81 


1.09 


+ 


0.19 


2.39 


± 


0.11 


10.6 + 0.7 


0.01 


± 


0.03 


1.13 


+ 


0.15 


189319 


-1.11 


+ 


0.13 


3978 


-0.29 


4.79 


+ 


2.03 


1.17 


+ 


0.17 


1.10 


± 


0.11 


48.3+5.1 


0.07 


+ 


0.04 


5.34 


+ 


0.89 


190608 


1.61 


± 


0.08 


4741 


0.05 


2.77 


+ 


0.94 


1.48 


± 


0.15 


2.82 


± 


0.08 


7.51+0.44 


-0.04 


+ 


0.03 


1.41 


+ 


0.13 


198232 


-0.45 


+ 


0.20 


4923 


0.03 


0.36 


+ 


0.07 


3.13 


± 


0.18 


2.47 


± 


0.07 


16.4+1.6 


-0.02 


± 


0.03 


1.30 


+ 


0.25 


198431 


1.45 


± 


0.17 


4641 


-0.12 


6.08 


± 


2.53 


1.17 


± 


0.14 


2.61 


± 


0.10 


8.53 + 0.76 


-0.03 


± 


0.03 


1.04 


+ 


0.17 


199665 


1.19 


± 


0.11 


5089 


0.05 


0.94 


+ 


0.18 


2.13 


+ 


0.13 


2.93 


± 


0.08 


8.00 ± 0.63 


-0.05 


+ 


0.03 


1.02 


+ 


0.13 


217428 


-0.24 


± 


0.33 


5285 


0.03 


0.36 


± 


0.08 


2.99 


+ 


0.21 


2.75 


± 


0.14 


11.6 + 1.9 


0.04 


± 


0.03 


0.55 


+ 


0.17 


218527 


0.75 


± 


0.24 


5084 


0.03 


0.79 


± 


0.19 


2.32 


± 


0.18 


2.81 


± 


0.06 


9.57 ± 0.85 


-0.08 


± 


0.02 


1.04 


± 


0.21 


219615 


0.68 


+ 


0.08 


4885 


-0.51 


5.46 


+ 


2.50 


1.03 


+ 


0.17 


2.43 


± 


0.06 


9.92 ± 0.42 


-0.01 


+ 


0.03 


2.30 


+ 


0.18 


224533 


0.70 


+ 


0.12 


5062 


0.00 


0.80 


± 


0.18 


2.31 


+ 


0.16 


2.80 


± 


0.04 


9.65 ± 0.54 


-0.05 


+ 


0.02 


1.31 


+ 


0.15 



